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Abstract 

^*\^ We study a single risky financial asset model subject to price impact 

r, and transaction cost over an infinite horizon. An investor needs to execute 

. a long position in the asset affecting the price of the asset and possibly 

C incurring in fixed transaction cost. The objective is to maximize the dis- 

C^-H counted revenue obtained by this transaction. This problem is formulated 

JL^ first as an impulse control problem and we characterize the value function 

I I using the viscosity solutions framework. We also analyze the case where 

there is no transaction cost and how this formulation relates with a singu- 
1/^ lar control problem. A viscosity solution characterization is provided in 

^ this case as well. We also establish a connection between both formula- 

tions with zero fixed transaction cost. Numerical examples with different 
types of price impact conclude the discussion. 



. Keywords: Price impact, impulse control, singular control, dynamic pro- 

r gramming, viscosity solutions 

o 



1 Introduction 

An important problem for stock traders is to unwind large block orders of shares. 



> 

X ^ , ^^ 

H Market microstructure literature has shown (e.g. Chan and Lakonishok 1995 

C^ 



Holthausen et al. 1990 ), both theoretically and empirically, that large trades 
move the price of risky securities either for informational or liquidity reasons. 
Several papers addressed this issue and formulated a hedging and arbitrage 
pricing theory for large investors under competitive markets. For example, in 



Cvitanic and Ma 1996 a forward-backward SDE is defined, with the price 
process being the forward component and the wealth process of the investor's 
portfolio being the backward component. In both cases, the drift and volatility 



'Department of Industrial Engineering and Operations Research, University of California, 
Berkeley CA 94720-1777. Email address: mjuncaOberkeley .edu 



M. JUNCA 



coefficients depend upon the price of the stocks, the wealth of the portfolio and 
the portfolio itself. Frey 1998 describes the discounted stock price using a 
reaction function that depends on the position of the large trader. In 



Bank 



and Baum 2004 Qetin et al. 2004 the authors, independently, described the 



price impact by assuming a given family of continuous semi-martingales indexed 
by the number of shares he ld ([Bank and Baum 2004] ) and by the number of 
shares traded ( [getin et aH|2004| ) 



The optimal execution problem has been studied in Bertsimas and Lo 1998 



Almgren and Chriss 2000 in a discrete-time framework. In both cases the 
dynamics of the price processes are arithmetic random walks affected by the 
trading strategy. In Bertsimas and Lo 1998 the impact is proportional to 



the amount of shares traded. In Almgren and Chriss 2000 the change in 



the price is twofold, a temporary impact caused by temporary imbalances in 
supply/demand dynamics and a permanent impact in the equilibrium or un- 
perturbed price process due to the trading itself. Also, this work takes into 
account the variance of the strategy with a mean-variance optimization proce- 



dure. Later on, nonlinear price impact functions were introduced in Almgren 



2003 . These ideas were adopted by more recent works under a continuous time 



framework. Schied et al. 2010 propose the problem within a regular control 



setting. The authors consider expected-utility maximization for CARA utility 
functions, that is, for exponential utility functions. The dynamics of the price 
and the market impact function are fairly general. Schied and Schoneborn 



2009 is the only reference that considers an infinite horizon model based on the 



original model in Almgren and Chriss 2000 . Finally, He and Mamaysky 2005 



consider only permanent price impact but they allow continuous and discrete 
trading (singular control setting) with a geometric Brownian motion as price 
process. 

On the other hand, it is also well established that transaction costs in asset 
markets are an important factor in determining the trading behavior of mar- 
ket participants. Typically, two types of transaction costs are considered in 
the context of optimal consumption and portfolio optimization: proportional 
transaction costs ( [Davis and Norniaii 1990 0ksendal and Sulem 2002] ) us- 
ing singular type controls and fixed transaction costs ( [Korn 1998 0ksendal 
and Sulem 2002] ) using impulse type controls. The market impact effect can 
be significantly reduced by splitting the order into smaller orders but this will 
increase the transaction cost effect. Thus, the question is to find optimal times 
and allocations for each individual placement such that the expected revenue af- 



ter trading is maximized. Ly Vath et al. 2007 include both permanent market 



price impact and fixed transaction cost and assume that the unperturbed price 
process is a geometric Brownian motion process. This reference only accepts 
discrete trading (impulse control setting) and uses the theory of (discontinuous) 
viscosity solutions to characterize the value function. Finally, [Subramanian 
and Jarrow 2001[ propose a slightly different model which does not include 
any fixed transaction cost but includes an execution lag associated with size of 
the discrete trades. It is important to remark that all papers referenced above 
assume a terminal date at which the investor must liquidate her position. 



OPTIMAL EXECUTION WITH PRICE IMPACT 



In this paper we study an infinite horizon price impact model that includes 



fixed transaction cost under the setting of impulse control, similar to Ly Vath 



et al. , 2007 . We describe a general underlying price process and a general per- 
manent market impact. With help of some classic results for optimal stopping 
problems and the discontinuous viscosity solutions theory for nonlinear partial 



differential equations, developed in references such as Crandall et al. 1992 Ishii 



and Lions 1990 Ishii 1993 Fleming and Soner 2006 , we obtain a fully charac 



terization of the value function when the price process satisfies some technical 
condition. Most of the processes used in financial studies satisfy this condition. 
This characterization is not complete when the fixed transaction cost is zero. 
By analyzing the HJB equation obtained in the previous case, we formulate a 
singular control model to include this case. For this new formulation we are 
able to complete the characterization. We are able to show that both formula- 
tions agree in the value function even though the formulations are completely 
different, when we choose the appropriate market impact functions. Finally, we 
describe the value function and the optimal strategy explicitly for an important 
special case. 

The structure of the paper is as follows: The general impulse control model, 
growth condition and boundary properties of the value function which are useful 
for the characterization of the function are exposed in Section 2. This section 
characterizes the value function of the problem as a viscosity solution of the 
Hamilton- Jacobi-Bellman equation and shows uniqueness when the fixed trans- 
action cost is strictly positive and the price process satisfies certain conditions. 
Section 3 proposes a singular control model to tackle the case when the trans- 
action cost is zero. Here a viscosity solution characterization and uniqueness 
result are proved as well under the same conditions. We present a special case 
where the value function of the impulse control model coincides with the value 
function of the singular control model and obtain the value function explicitly 
for this case. Section 4 shows that, in fact, these two formulation produce the 
same value function even though they consider different types of control. Sec- 
tion 5 presents numerical results for different underlying price processes with 
both formulations. Finally, we state some conclusions and future work. 



2 Impulse control model 

Let (rj, J^, {J-t)o<t<oa,P) be a probability space which satisfies the usual condi- 
tions and Bt be a one-dimensional Brownian motion adapted to the filtration. 
We consider a continuous time process adapted to the filtration denoting the 
price of a risky asset Pt- The unperturbed price dynamics, when the investor 
makes no action, are given by: 

dPs = ^i{Ps)ds + a{Ps)dBs, (2.1) 

where ^ and a satisfy regular conditions such that there is a unique strong 
solution of this SDE (i.e. Lipschitz continuity). We are mainly interested in 
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price processes that are always non-negative, thus we assume that P is absorbed 
as soon as it reaches and that the initial price p is non-negative. 

Now, the number of shares in the asset held by the investor at time t is 
denoted by Xt and it is up to the investor to decide how to unwind the shares. 
Different models and formulations will define the admissible strategies for the 
investor. At the beginning the investor has x > number of shares and we only 
allow strategies such that Xt > for all t > 0. Since the investor's interest 
is to execute the position, we don't allow to buy shares, that is Xt is a non- 
increasing process. Hence, we can see that K+ x R+ = O (with interior O) 
is the state space of the problem. The goal of the investor is to maximize the 
expected discounted profit obtained by selling the shares. Given y = {x,p) G O 
we define V{y), the value function, as such maximum (or supremum), taken 
over all admissible trading strategies such that (Xo_, Pq-) = Yo- ~ V- We call 
/3 > the discount factor and fc > the transaction cost. Note that we can 
always do nothing, in which case the expected revenue is 0. Therefore ^ > 
for all y. 

In this formulation we assume that the investor can only trade discretely over 
the time horizon. This is modeled with the impulse control v — {Tn, Cn)i<n<M, 
where the random variable M < oo is the number of trades, (r„) are stopping 
times with respect to the filtration (Tt) such that < Ti < • • • < t„ < • • • < 
tm < oo that represent the times of the investor's trades, and (Cn) are real- 
valued J>^ -measurable random variables that represent the number of shares 
sold at the intervention times. Note that any control policy v fully determines 
M. Given any strategy ly, the dynamics of X are given by 

Xs = Xr^, for Tn<S< T„+i, (2.2) 

^T„ + i = Xt^ — Cn+1- (2-3) 

We consider price impact functions such that the price goes down when the 
investor sells shares. Also, the greater the volume of the trade, the grater the 
impact in the price process. Thus, we let a{(,p) be the post-trade price when 
the investor trades ( shares of the asset at a pre-trade price of p. We assume that 
a is smooth, non-increasing in (, and non-decreasing in p such that a{0,p) — p 
for all p. These conditions imply that a{(,p) < p for C ^ 0- Furthermore, we 
will also assume that for all Cij C27P G 1^+ 

a(Ci,a(C2,p))-a(Ci + C2,p). (2.4) 

This assumption says that the impact in the price of trading twice at the same 
moment in time is the same as trading the total number of shares once. This 
assumption will prevent any price manipulation from the investor. Two possible 
choices for a are: 

aiiCp) =P- K 

a2{C,p) =pe^^'^ 

where A > 0. A linear impact like ai has the drawback that the post-trade price 
can be negative. Given a price impact a and an admissible strategy v, the price 
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dynamics are given by: 

dPs = fJ.{Ps)ds + a{Ps)dBs, for t^ < s < r„+i, 
Pr„ =a(C„,Pr„J. 

When T„ = T„_|_i, then we apply the impact twice, therefore 



(2.5) 
(2.6) 



Pr„ = Pt„ + i = a{Cn+l,a{Cn,Pr„^))- 

If more that two actions are taken at the same time, we apply the impact 
accordingly. Now, given y — {x,p) G O the value function V has the form: 



V{y) = supE 



M 



-/Jr„ 



(CnPr„ - k) 



(2.7) 



As usual, we assume that e ^'^ = 0on{r = o3}. 

2.1 Hamilton- Jacobi-Bellman equation 

In order to characterize the value function we will use the dynamic programming 
approach. This principle has been proved for several frameworks and types of 
control. Some of the references that prove it in a fairly general context are 



Ishikawa 2004 Ma and Yong 1999 . We have that the following Dynamic 



Programming Principle (DPP) holds: For all y = {x,p) £ O we have 



V{y) = supE 



E' 



.-I3t„ 



{CnPr^-k)+e-^^V{Yr) 



(2.8) 



where r is any stopping time. Let's define the impulse transaction function as 

r(2/,C) = (a^-C,a(C,p)) 

for all y G O and C G K. This corresponds to the change in the state variables 
when a trade of ( shares has taken place. We define the intervention operator 
as 

M^piy)^ sup (p(r(y,C)) + Ca(C,p)-fc, 

0<C<2: 

for any measurable function ip. Also, let's define the infinitesimal generator 
operator associated with the price process when no trading is done, that is 



dp 



A^ = ^i{p)^ + ^<j{pf 



ay 

dp'^ ' 



for any function (p G C'^{0). The HJB equation that follows from the DPP is 
then ( |0ksendal and Sulem[|2005| ) 



min {(3ip — Aip, tp — Ai(p} = in C 



(2.9) 
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We call the continuation region to 

and the trade region to 

T={yeO ■.Mip-ip = 0}. 



2.2 Growth Condition 

We will define a particular optimal stopping problem and use some of the re- 
sults in [Dayanik and Karatzas^ ,2003] to establish an upper bound on the value 
function V and therefore a growth condition. Consider the case where there is 
no price impact, that is, a{(,p) = p for all ( > 0. We define 



VNiiv) =supE 



M 



J2e-^^-iCnPr„-k) 



(2.10) 



where Pg follows the unperturbed price process. It is clear that V < Vni- When 
there is no price impact, the investor would need to trade only one time. 



Proposition 2.1. For all y ^ O 

Vni{x,p) = U{x,p) := supE[e-'^^(2;P^ - k)+] 



(2.11) 



where the supremum is taken over all stopping times with respect to the filtration 

Proof. Since (r, cc) is an admissible strategy for any stopping time r, then U < 
V/v/- Now, let T„ the set of admissible strategies with at most n interventions. 
The proof will continue by induction in n to show that for all n 



sup E 



Y,e-^^HCrPr.-k) 



< U{y). 



(2.12) 



Clearly (2.121 is true for n — 1. Let u e T„. Note that xp — k < U{x,p), 
therefore, conditioning on J-V^ we have 



E 



J2e-^^'{C^Pr.-k) 



.i=l 



=E[E[[e-^-HCiPr,~k)\^r,]] 



E 



E 



-/3ti 



E' 



'|3{r^-Tl) 



{QPr. - k) 



Tr 



<E[E[e"'^^^C/(Cl,P.J|-F,J]^ 
E [E [ e^'^^i C/(a; - Ci , Pri ) I Tr 
<E[e-^-iC/(a;,P,J] 
<V{x,p), 
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where the last inequahty fohows from the fact that the process e~^'^U{x,Ps) 
is a supermar tingale ([0ksendal and Reikvam 1998| ). This proves ( |2.12[ ). By 
Lemma 7.1 in 0ksendal and Sulem 2005 , the left hand side of (2.12) converges 
to Vni as n — )• oo and the proof is complete. D 



oo and the proof is complete. 
From the previous result we have the bound 

< V{x,p) < U{x,p) = supE[e"'^^(xP^ - fc)], 



(2.13) 



where the supremum is taken over all stopping times with respect to the filtra- 
tion (J-t). Following section 5 in [Dayanik and Karatzas 2003 , let ip and 4> be 
the unique, up to multiplication by a positive constant, strictly increasing and 
strictly decreasing (respectively) solutions of the ordinary differential equation 
Au — (3u and such that < '0(0+) and ^{p) — >■ oo as p — >■ oo. For any x > 0, 
let 

4 = lim 1^E_^. (2.14) 

Then U is finite in O if and only if i^x is finite for all a; > 0. Furthermore, when 
U is finite wc also have that for some C > 



and 



U{x,p) < Cx^pip) 



lin, ^^ ^ 4 



p-J-oo 1p(j)) 



(2.15) 
(2.16) 



We will assume that U is finite. 



2.3 Boundary Condition 

Since the investor is not allowed to purchase shares of the asset we have that 
V{0,p) = for all p > 0. Also, the price process gets absorbed at 0, there- 
fore y = on do. If we assume that U is finite then by ( 2.15[ ) we have that 
V{x,p) — > as x — >■ for all p > 0, that is, V is continuous on {x ~ 0}. Now 
we distinguish two cases: 



1. is an absorbing boundary for the price process P. This means that for 
any p > 0, P{Pt = for some t > 0\Po = p) > 0. A simple example is the 
arithmetic Brownian motion. Since the process is stopped at 0, we must 
have that for all a; > 

U{x,0)^0. 

Also, [Dayanik and Karatzas[ 2003| shows that in this case U is continuous 
at {p = 0} whenever U is finite. Therefore the boundary conditions for 
the value function V are 



V = 0ondO and lim V{y') 



for all y £ dO. 



(2.17) 
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2. is a natural boundary for the price process P. This means that for any 
p > 0, P(Pt = for some t > 0\Po = p) = 0. For example the geometric 
Brownian motion. In this case we can have different situations in V{x,p) 
as p goes to depending on the price process. In particular, we can have 
the situation where V is discontinuous on the set {p = 0}. 

2.4 Viscosity solution 

We now are going to prove that the value function is a viscosity solution of the 



HJB equation (2.91 and find the appropriate conditions that make this value 



function unique. The appropriate notion of solution of the HJB equation (2.9 1 
is the notion of discontinuous viscosity solution since we cannot know a priori 
if the value function is continuous in O. We must first state some definitions. 

Definition 2.2. Let W be an extended real-valued function on some open set 

(i) The upper semi-continuous envelope of W is 

W*{x) = lim sup W{x'), Vx G V. 

(a) The lower semi-continuous envelope of W is 

WJx) = lim inf W(x'), Va; € V. 

riO |x'-x|<r 

x'e-D 

Note that W* is the smallest upper semi-continuous function which is greater 
than or equal to W, and similarly for W^,. Now we define discontinuous viscosity 
solutions. 

Definition 2.3. Given an equation of the form 

mhi{F{x,(p{x),D(p{x),D^ip{x)),ip-Mip} =0 inV, (2.18) 

a locally bounded function W on V is a: 



(i) Viscosity subsolution of (2.18) in T) if for each Lp e C^{'D), 

min {F{xo, W{xo), D^{xo), D^ip{xo)), W*{xo) - MW*ixo)} < 

at every xq (z T) which is a maxiniizer of W* — Lp on T) with W* {x^)) 
Lp{xo). 



(a) Viscosity supersolution of (2.18) in V if for each ip G C {V), 

min {F{xo, W[xo),D^{xo),D\{xo)).W,{xo) - MW,{xo)) > 
at every xq Cz V which is a minimizer ofW^,—p onV with W^,{xq) — ip{xo). 
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(Hi) Viscosity solution of (2.18) in T) if it is both a viscosity suhsolution and 
a viscosity supersolution of (2.18) in V. 



We are now ready for the following theorem: 



Theorem 2.4. The value function V defined by (2.7) is a viscosity solution of 



(2.91 mO. 



Proof. By the bounds given in the section |2.2[ it is clear that V is locally 
bounded. Now we show the viscosity solution property. 

Subsolution property: Let yo ^ O and </? G (7^(0) such that j/o is a maximizer 
of V* — (p on O with V*(yo) ~ fiua)- Now suppose that there exists 9 > and 
(5 > such that 

- /3^{y) + A^{y) < ~9 (2.19) 

for all y e O such that \y — yo\ < S. Let (?/„) be a sequence in O such that 
Vn ~> 2/0 and 

lim V{y^)^V*{yo). 



By the dynamic programming principle (2.8 1, for all n > 1 there exist an ad 



missible control Vn = i'''mTCm)m such that for any stopping time t we have 
that 



^(2/n) < E 



J2 e-^^-(CF;v -k) + e-^-i/(i;") 



1 



(2.20) 



where YJ^ is the process controlled by z/„ for s > 0. Now consider the stopping 
time 

Tn^mi{s>0:\Yr~yo\>S}ATl\ 



where rf is the first intervation time of the impulse control f„. By (2.20) we 
have that 

v{yn) < E [e-^^-v{Y^\j^T,.<rn] +E [e"^^" (cp."^ - k + y(r;f )) iiT,.=rn 



< E [e-^^"F(r£^_)i{T„<.r}] +E[e-^^"A^y(r;i_)i{T„=.n] + 1 

< E [e-^^-V{Yil^)] + - 



(2.21) 
(2.22) 



1 

n 



Now, by Dynkin's formula and (|2.19|) we have 



E[e-'^^"(^(r^^_)]=(^(y„)+E 



e-^^-MYn+AipiYr))ds 



<^{y^)--{l-E[e-"^-]). 
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Since V < V* < ip and Tn < rf, by (2.221 



p n 

for all n. Letting n go to infinity we have that 

lim E[e"'3^"] = 1, 

n— >-oo 

which implies that 

hm P[ti" = 0] = 1. 

n— j-oo 



Combining the above with (2.21 ) when we let n — )- oo we get 



V*{yo)< sup MV{y'). 
\y'-yo\<S 

Since this is true for all S small enough, then sending (5 to we have 

V*iyo)<iMVnyo). 

If we show that {MV)* < MV* , then we would have proved that if —(3(p{yQ) + 
Aip{yo) < 0, then MV*{yo) - V*{yo) > and therefore 

min{/3(p(yo) - Aip{yo),V*iyo)-MV*iyo)} < 0. 

Appendix |A] contains the proof of this last fact. 

Supersolution property: Let y^ G O and (p G C^{0) such that j/o is a 
minimizer of T4 — (/j on O with Kt(i/o) = "i^lyo)- By definition of V and A^F we 
have that M.V < F on O and therefore {^AV).^, < V^,. Let (y„) be a sequence 
in O such that y„ — ?► yo and 

lim V'(2/„) = K(yo)- 

Now, since K^ < 1/ is lower semi-continuous and F is continuous we have 

MV,{yo)^ sup V4T{yo,0) + CaiC,po)-k 

0<C<xo 

< sup liminfF(r(y„, C)) + Ca(C,Pn) - fc 
<liminf sup F(r(y„, C)) + Ca(C,Pn) " ^ 

n-foo o<C<a;„ 

< lim MV{yn) 

n—>-oo 

= {MV),{yo). 

Hence MV^,{yo) < {MV)^,{yo) < K(yo)- Now suppose that there exists 9 > 
and (5 > such that 

/3(p(y) - Acpiy) < -0 (2.23) 
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for all y G O such that \y — yo\ < S. Fix n large enough such that |?/„ — yo| < ^ 
and consider the process YJ^ for s > with no intervention such that Fq" = 2/„. 
Let 

T„=inf{s>0:|n"-2/o|>'^}- 



Now, by Dynkin's formula and (2.23) we have 



El 



-0T„ 



^(y^J]=^(y„)+E 



'Ps 



{~MYn+A^iYr))ds 







><p(y„) + -(l-E[e-^^"]). 
On the other hand, if <V^, <V and using the dynamic programming principle 



(2.81 we have 



Notice that 77 := lini E[e^'^"^"] < 1 since T„ > a.s by a.s continuity of the 

n— >-oo 

processes YJ^, then by the above two inequalities and taking 7i — >■ 00, we have 
that 

K(yo) > (p{yo) + -^{l-r])> (p{ijo) 

contradicting the fact that Kt(yo) = ¥'(yo)- This establishes the supersolution 
property. D 



2.5 Uniqueness 



Let -0 be defined as before and let's assume that the function U defined in (2.131 



is finite. Also assume that the transaction cost fc > 0. Then, we want to prove 



that V is the unique viscosity solution of the equation (2.9 1 that is bounded by 



U. We will need an additional assumption about the function -0: For all a; > 



hm = £^ 

p->-oo 1p[p) 



0. 



(2.24) 



Following the ideas in [Crandall et al. [1992 Ishii 1993 let u be an upper semi- 
continuous (use) viscosity subsolution of the HJB equation (2.9) and u be a 
lower semi-continuous (Isc) viscosity supersolution of the same equation in O, 
such that they are bounded by U and 



limsupw(y') < liniinf w(y') for all y G dO. 



v'->-v 



(2.25) 



Define 



Vm(x,p) = v{x,p) H X i){p) 

m 
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for all m > 1. Then Vm is still Isc and clearly (3vm — Avm > by definition of 
•0. Now, 

Mvjn{x,p)= sup v{x-C,a{(,p))-{ {x~Cftlj{a{(,p)) + Ca{(:;,p)-k 



0<C<2; 



1 



< sup v{x-C,a{C,p)) + Ca{C,p)-k+ sup — (x - C)V("(C,p)) 

0<C<a; 0<C<2: "^ 

= A^u(a;,p) H x'^^{p) 

m 

< v{x,p) H X^tpip) = V,n{x,p). 

m 



V and equations (2.151 and (2.24) we get 



Therefore Vm is supersolution of (2.9). Now, by the growth condition of u and 

(2.26) 



\y\ 



lim {u-Vjn){y) = -oo. 



We will show now that 



u < V in O. 



(2.27) 



It is sufficient to show that sup(?i — w.,„) < for all to > 1 since the result 

yeo 
is obtained by letting to, — > cx). Suppose that there exists m > 1 such that 

?7 — sup(m — Vm) > 0. Since u — Vm is use, by (2.26) and (2.25) there exist 

yeo 
yo (z O such that rj = [u — Vm){yo)- Let yo = {xotPo) be the one with minimum 
norm over all possible maximizers of u — «,„. For J > 1, define 

My,y') = 2^y-y'\'^ + \y-yo\'^, 
'^i{y, y') = u{y) - vm{y') - My^ y')- 



Let 



sup ^i{y,y')^^i{y,,y[) 

\y\,\y'\<\yo\ 

'|4 



Clearly -q^ > 77. Then, this inequality reads ^\yi - y-p + jy,, - j/ol < "(2/») - 
Vm{y'i) - (u - Vm){yn)- Since |yi|,|2/-| < |yo| and u and -«„ are bounded 
above in that region, this implies that 2/^,2/^ — ^ yo and ||yi — J/^l'* — >■ (along 
a subsequence) as i — 7> 00. We also find that rji — >■ 77, u{yi) — Vm{y'i) — >■ t? and 
u{yi) -^ u{yo),Vrn{yi) -^ viyo)- By theorem 3.2 in Crandall et al.[ 1992] , for all 
i > 1, there exist symmetric matrices Mi and M- such that {^^{yi,y[),Mi) = 

{d,,Mi) e P'+u{yi), {-^{y^,y^,Ml) = {d'^MD e p.-^^(yo'and 



M, 
M 



,) <D^My^,y'i) + ]iD'My^,y'i}y 



Since u is a subsolution of (2.9) and Vm is a supersolution, we have 

mm{f3u{yi) - n{pi)di^2 - -cr{pif Mi ^22, u{yi) - Mu{yi)} < 0, 



OPTIMAL EXECUTION WITH PRICE IMPACT 



13 



and 



1 



{/3v„M) - /^(pD^b - o'^(P^)'K22,«m(yD - Mv^iy',)} > 0. 



Now, if we show that for infinitely many Vs we have that 

1 



and since it is always true that 



-a{p,YM,^22 < 0, 



<^{p'ifM' > 0, 



(2.28) 



we have that u < Vm by following the classical comparison proof in Crandall 



et al. 1992 . Suppose then, that there exists io such that (2.28) is not true for 



all i > io, then for i > io 

u{yi)~Mu{y,) < 0. 

Since Vm is a supersolution, we must have that 

v^{y[)-Mv,n{y'i)>0. 

Since u is use, there exist Q such that A4u{yi) = u{xi — Q, ctiCii Pi)) +Ci'^iCi^ Pi) ~ 
k. Then 

u{yi) < u{xi - Ct,a{Ci,Pi)) +Cia{Q,Pi) - k. 

Extracting a subsequence if necessary, we assume that Q ^ C,q as i —^ oo. 
First, consider Co = Oi then by taking limsup in the inequality above we get 
u(jJo) < u(yo) — k. This is a contradiction since fc > 0. Now assume that Co =/= 0. 
From the above inequalities we have that 

u{yi)~Um{yl) < u{x^~C^,a{Ci,Pi))+C^(^{Ci,Pi)-Vm{x'^~C^,a{Cip'i))-C^(^iC^,P'i), 

for any < Cl < Pi- Since p'^ — >■ Poj let Ct ~^ Co ^-^d taking limsup in the above 
inequality we get that 

V < {u-Vrn){xQ -Co,a(Co,Po))- 

This is a contradiction since j/o was chosen with minimum norm among maxi- 
mizers of u — Vm and Co > 0. Therefore (2.28) must hold for infinitely many I's 



and (2.27) holds. As usual continuity in O and uniqueness of V follow from the 



fact that T^ is a viscosity solution of (2.9) 



We have just proved the following theorem: 



Theorem 2.5. Assume condition (2.24) and that the transaction cost fc > 0. // 



W is a viscosity solution of equation (2.9) that is bounded by U and satisfies the 



same boundary conditions as V, then W = V. Furthermore, V is continuous in 
O. 



Remark 2.6. Condition (2.24) is satisfied by ltd processes like Brownian Mo- 



tion, Geometric Brownian Motion, Ornstein- Uhlenbeck and Cox-Ingersoll-Ross. 
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3 No transaction cost 

From the proof of the above uniqueness resuh, we can see that the resuh de- 
pends on the fact that fc > 0. Let's start by pointing out that in this case the 
intervention operator becomes 

Mifiy)^ sup ip{Tiy,C)) + CaiC,p)>viT{y,0))^ipiy), (3.1) 

for any measurable function (p. This imphes in particular that any measurable 



function is a viscosity subsolution of (2.9). On the other hand, V > MV for 
the value function. Then we have that 

V > MV > V. 

Assume now that V G C^{0). Since C = is a maximum for ( i-^ V{T{y, ()) + 
Ca{(,p), then for aU j; e O: 



C=o 



da dV dV da 

da dV dV 

Recall that a is non-increasing in ('j so we define 

l{p) = ~{0.p). (3.2) 

for all p > 0. Hence, we get the following condition for V: 

dV dV 

-i{p)^iy)-^{y)+P<0- (3.3) 

dp dx 

This suggests that if we assume no fixed transaction cost we should look at a 
different HJB equation, that is 

mm<;/3(^-Av3,7(p)^ + ^-pl =0. (3.4) 



On the other hand, condition (2.4) implies that it is always better to split the 



orders into smaller orders. Indeed, given (C,p) G O and < C' < C 

C'a(C',p) + (C - C')«(C,P) = Ca(C,p) + (C - C')("(C',p) - «(C,p)) > C«(C,p), 

since a is non-increasing in Q. 
3.1 Singular control 



In fact, the equation (3.4) is the associated equation of the following control 
problem ([Oksendal and Sulem 2005|): In this case our admissible controls are 



of the singular type, that is 

dXt = -diu 
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where ^o = 0, <^ is an adapted, continuous non-decreasing and non-negative 
process. The price process in this case foUows the dynamics 

dPt = fiiPt-)dt + o{Pt-)dBt - -l{Pt-)diu 



where 7 (see (3.2 )) is a non-negative smooth function that accounts for the price 
impact. In order to guarantee the existence and uniqueness of the process Pt, 



we need to also assume that 7 is a Lipschitz function ( Protter 2004| ). Now, 
the form of the value function Vq changes to 



Vo{y) 



supE 



.-pt 



Ptd^t 



for all y € C In this case the appropriate form of the DPP is 



K)(y) = supE 



.-/9s 



Psd^s 



.-Pr 



Vo{Yr) 



(3.5) 



(3.6) 



for any stopping time r. As before, we can define the continuation region as 

dip dip 



C^{y£0: 7(p) 
and the trade region as 



dp 

dip 



dx 



dp 



dp dx 



p>0} 



-p = 0}. 



Typically, singular controls are allowed to be cadlag instead of continuous. We 
decide to restrict our controls for two reasons: (1) Under the absence of fixed 
transaction cost, the investor will divide the orders into very small pieces as 
shown above. (2) When the singular control is discontinuous the stochastic 
integral may not be properly defined (see Protter 2004|). 



3.2 Viscosity solution 

Although we only consider continuous strategies, the value function is still a 



viscosity solution of equation (3.4) (which definition is similar to 2.3 1. 
Theorem 3.1. The value Junction Vq defined by (3.51 is a viscosity solution of 



(3.41 mO. 



Proof. Since we can approach finite variation functions by simple functions, by 
proposition |2 . 1 1 we have that 

Vo < U. (3.7) 

Therefore, Vq is locally bounded. 

Subsolution property: Let yo ^ O and ip G C^(0) such that yQ is a maximizer 
of Vq — ip on O with yQ*(yo) ~ fiva)- Now suppose that there exists k > and 
6 > such that 

- Pviy) + Ap(y) < -K and p - -i{p) — {y) - —{y) < -n (3.8) 
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for all y G O such that \y — yo\ < S. Let (y„) be a sequence in O such that 
J/n -^ 2/0 and 

lim Vb(2/„) = VQ{yo). 



Given any stopping time r, by (3.6), for all n > 1 there exists an admissible 
control ^" such that 



1^0 (yn) < E 



Jo 



1 

n 



where Y" is the process controlled by ^" for s > starting at y„. Since Vq < 
Vo — fy using Dynkin's formula for semimartingales ( [Protter 2004 ) we have 

that 



Vo(yn) < E 

-E 



Jo 



ViVn 



-'' hiPll^iYI' 



-E 
dx 



e-P^~l3^{YP) + A^{Y:))ds 



(YD dC 



1 

n 



Consider again the stopping time 

T,,=mi{s>0:\Yr-yo\>S}, 



then by (3.8 1 



v^o(yn) < -'^E 



e-^^ids + dC) 



+ 'fiiVn) + - 



Taking n — >■ c» we obtain a contradiction since the integral inside the expecta- 
tion is bounded away from for any admissible control ^ by the a.s continuity 



of the process Yp. Hence at least one of the inequalities in (3.8) is not possible 
and this establishes the subsolution property. 

Supersolution property: Let yo E O and ip E C^{0) such that j/o is a 
minimizer of Vo* — ip on O with Vo*(yo) = 'fiiya)- Let (?/„) be a sequence in O 
such that yn -^ j/o and 

lim V"o(y„) = V'o^yo)- 



First, suppose that there exists 9 > and S > such that 

I3^{y) - A^{y) < -9 



(3.9) 



for all y e O such that \y — yo\ < S- Fix n large enough such that \yn — Vol < S 
and consider the process Y" for s > with no intervation, i.e. ^ = 0, such that 
Fo" = y„. Let 

T,,=mi{s>0:\Y,^-yo\>S}. 
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Now, by Dynkin's formula for scmimartingales and (3.91 we have 

-E 



e-P^{-(3ipiYr) + A^{Yn)ds 



^^-^M7(pr)|(>^.") + g(n"))rf6 



= ifiyn) + E 
> ifiy„) - eE 



e-^^{^/3ipiYr) + A^{Yn)ds 



.-Ps 



ds 



As before, from here we can draw a contradiction with Vb*(yo) = 'fivo) by the 
a.s. continuity if the process Y^ . Now, take /i > and consider the process Yt 



with control process d^t — ji^[o.h]{t)dt and Fq = y for given y £ O. Using (3.6 1 
we can show that 



Vo{y) > E 



>E 



e-P'Psd^s+e-P'^ViYh) 



/3hT 



■ nh 

Jo 



-Ph 



viYh] 



= E ' 
By Dynkin's formula again, 

E[e-^V(n)]=¥'(y)+E 



-/3s 



P,ds 



-Ph 



^{Yh) 



-^'{-l3^{Y,) + A^{Y,))ds 







E 



''(iiPsf^Ays) + 'liYMdi. 







dp 



dx 



^(zj)+E 
1 



-'''{~P^{Ys) + A<p{Ys))ds 



-E 



r.-(.«)|«,. !«,).. 



Letting /i — >• 0, we have 



1 o 



Therefore, for all n > 1 we have 



Vb(yn) > (fiiyn) +p,i - i{pn)-Q-{yn) - -Q^iy^)- 
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Since 7 is continuous, letting n — >■ 00 we get 

as desired. This establishes the supersolution property. 



D 



3.3 Uniqueness 

Recall that with the impulse formulation we do not have uniqueness in absence 
of transaction cost. This is not the case with the singular control formulation. 



Theorem 3.2. Assume that (2.24 1 is satisfied. If W is a viscosity solution of 



equation (3.4 1 that is bounded by U and satisfies the same boundary conditions 



as Vq, then W = Vq. Furthermore, Vq is continuous in O . 

Proof. The proof follows the same strategy as in the impulse control case. Let 
u be an upper semi-continuous (use) viscosity subsolution of the HJB equation 
(3.4) and w be a lower semi-continuous (Isc) viscosity supersolution of the same 
equation in O, such that they are bounded by U and condition (2.25) holds. 
Define 

Vrnix,p) = (1 ) V{x,p) + — {C{x + lf^j{p) + l) 

\ 771/ 777 



for all 7n > 1 and C as in (2.15). Recall that 7 is non-negative and ip is an 
m, 



increasing function, then (2.15) implies that 
dv. 



-P 



dx 



lip)- 



dp 



> - 



P 



l--)p+|--C(x-Kl)V(p)+7(p)|--C(x+l)V(p) 
rri J ax m opm 



> 



-p + -2C{x + l)7/'(p) + j{p)-C{x + 1) Vb) 

777 m, 777 

12 1 

— P+ —p{x + l)+'j{p) — C{x + lfilj'{p) 
mm m 



1 

>-P- 

777 



_ /3 



Also(/3/-A)(i) 

a strict supersolution of (|3.4|) in O. Following the same lines and definitions as 



> 0, where / is the identity operator. Therefore Vm is 



in the previous proof we have 



mm{l3u{yi) - p.{pi)d^,2 - -a-{pt)^Mi^22,-Pi + di,i + liPi)di,2} < 0, 



and 



mm{f3v,n{y',) - p{p'i)4^2 - ^'^{p'i)^K22, -P\ + ^b + 1{P'^¥^,2} > ^^. 

where bi = niin i ^, ^ f • Since p\ — >■ po and 7/0 G Cj ^i > for large enough i. 
We need to show now that for infinitely many i's we have that 



/37t(j/i) - p.{pi)di,2 - -o(jliYMi^22 < 0. 



(3.10) 
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Suppose then, that there exists «o such that (3.101 is not true for all i> ia, then 
for i > io 

-Pi + dis + 7{Pt)di,2 < 0. 

Since Vm is a supersolution, we must have that 

-p'^+dil+l{p'M^2>S^. 

Hence, 

Pi-P'z - {d^^ - d'i{) - {-f{pi)d,^2 ~ l{p'i)d'i,2) > ^^■ 

Since di,d'^ goes to as i goes to oo, we get the contradiction > S^ — 
min -^ m' m ( -^ ^- Therefore (3.10) must hold for infinitely many I's and the 
comparison result holds. Everything follows now as before. D 



3.4 Optimal strategy for a special case 

Previous sections characterize the value function of our problem in different 
formulations. We will calculate now the explicit solution of the value function 
and describe the optimal strategy in a particular case. Let us come back to the 
impulse control case. Since we are allowed to do multiple trades at the same 



time, we are going to explore this strategy. Assumption (2.4) guarantees that 



the price impact does not change by splitting the trades, but the profit obtained 
by doing so could be greater. Let's define the following function 



W{y) 



/ a{s,p)ds for y £ O. 
Jo 



(3.11) 



This is the best that we can do when we do many trades at the same time. It 
is clear that this is not attainable with any impulse control. Since a is non- 
increasing on X and positive, we have for all y G O 



xa{x,p) < / a{s,p)ds. 
Jo 



(3.12) 



Therefore, for all < C < 2; 

W{T{yX)) + CaiCp) = Ca{C,P) + 
<Ca(C,P) + 



x-c 



a{s,a{C,p))ds 



C 



a{s,p)ds— / a[s,p)ds 
^0 

= W{y) + Ca{C,p)- a{s,p)ds 
Jo 

< W{y), 



where the last inequality follows from (3.12). Hence MW <W and therefore 



MW = M^ by (3.1). On the other hand, the function W satisfies (3.3) with 
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equality. Indeed, by the condition (2.4) we have that for any Ci, C2 and p 
da , 



dc 



(Cl+C2,p) 



-Q-{Ci,a{C2,p))-^{C2,p), 



and taking C2 = we obtain 



9a, _, . da \9a,^ . , .da,^ , 



Now, since a is smooth we find 
dW 



dW 



-Q^iy)+p = ~iip) -Q-is,p)ds 



da 



dx 



a{s,p)ds + p 







"9a 



{s,p)ds — a{x,p) +p 



dC 
a{x,p) — a(0,p) — a{x,p) +p - 



0. 



If we had also that [3W — AW > 0, then W would solve both equations (2.9 1 
and (l33 and r=0. 



Now, Subramanian and Jarrow 2001 considers impact functions of the form 



a{x,p) = pc{x), where < c < 1 is nonincreasing. In our case, by condition 
(2.4), c must satisfy c(xi)c{x2) = c{xi + X2) and therefore we end up with the 
following price impact functions and W: 



a(x,p) = pe 
l{p) = Ap 
W{x,p) = ^(1 



■\x 



-\x\ 



(3.13) 
(3.14) 

(3.15) 



with A > 0. This function was proposed also in He and Mamaysky 2005 and 



Ly Vath et al. , 2007 . Let's consider this price impact function for the moment. 



In this case we have the following: 

Theorem 3.3. V^^W = V if and only if U{x,p) 



xp. 



Proof. If Vb = M^ then PW - AW > and therefore (3ip- Aip>0 for ip{p) = p. 
By the uniqueness result for optimal stopping problems (see Theorem 3.1 in 



0ksendal and Reikvam 1998 ) 



p = supE[e-'^'"P,-], 

r 

that is U{x,p) — xp. Suppose that 

U{x,p) — xsupE[e^^'^P^] = xp, 



for y e C This means that I3ip-Aip>0 for (j){p) = p. Therefore fiW - AW > 
and W satisfies the HJB equation (3.4) with T = O. Also, W satisfies the 
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growth condition and has the same boundary conditions as Vq by (2.13). By 



Theorem 3.2 we have that W = Vo- To prove the second equahty we will do 



induction in the number of trades. Note that the function C i— ^ ^e '» in [0 
attains its maximum at a; = minja::, j}. Then, 

sup E[e-^''^CiPr,^e-^'^^] < U{x,p) = xp < W{x,p). 
Now, let v e T„. Hence, 

On the other hand, by induction hypothesis we have 



E 



-/3n^g-0(r.-n)(>^p^^ 



= E 



-/3ri 



E 



^e-^(--^)C.P.. 



Tr 



<E[e-^^'E[Wix - Ci,e-^'^'Pr,-)\Tr,]] 



Combining both inequalities above we have 



E 



E' 



-/3r. 



c^K 



< E[e~^'''W{x,Pr,-)] < Wix,p). 



Again, by Lemma 7.1 in 0ksendal and Sulem 2005 , the left hand side converges 



to F as n — )■ cx). Clearly the other inequality holds and the proof is complete. D 
Example Consider the case where the price process is a geometric Brownian 



motion. This is the only process that is considered in the papers Subrama- 



nian and Jarrow 2001 He and Mamaysky 2005 Ly Vath et al. 2007 . The 



unperturbed price process is 

dPt = fiPtdt + crPtdBt, 

with cr > 0. It is easy to see that the value function U is finite if and only if 
/3 > /i. In this case the function i/j takes the form 

where v > 1, therefore condition ( |2.24[ ) holds. Now, the condition (2.13) reads 

< V{x,p) < U{x,p) = xp. 

This implies that Vq = V — W. We can see that in this case the value function 
W is not attainable with any impulse control, but we can approach it by trading 
smaller and smaller orders. We will show now how we can approach W with 
singular (in fact regular) controls. Let u > and consider the strategy d^t = udt, 
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that is, selling shares at a constant speed u until the investor executes the 
position. Then, 



and 



Pt = pexp{(^ - Am - -(T^)t + aBt} 



E 



e-P'Ptd^t 



= vE 



x/u 





x/u 



e -'^'Ptdi 



w / e 
'o 

up / e('^-^"''5)*di 

"'0 



pu 



Jp-Xu-P)x/u 



1 



fi — \u — j3 \ 

by using Fubini's theorem since the integrand is positive. Taking u ^ oo this 
expression converges to W . 



4 Connection between both formulation 



Theorem 3.3 shows that V — Vq for a special case, i.e., the value function of 
two different problems are the same. We are going to show that this is not a 
coincidence. Let us start with some notation: Given fc > and y ~ {x,p) (z O 
we denote: 

M 



X/W(y)=SUpE 



Y^e'^^-iCnPr^-k) 



.n=l 



and 



x(^-V(2/)= sup ^(r(y,c)) + C«(C,p)-fc- 

0<C<2; 



Lemma 4.1. For all y E O we have 

limy('=)(y) = F(°)(y). 

fc— ^0 

Proof. It is clear that V*^"-' is an upper bound. Let e > 0, then there is m > 
and v G T„i such that 



l/(")(2/)<E 
For any fc < — we have that 

1/(0) (y) <E 

<V'^^>iy) + 2e 



Y.e-^^'Qa{Q,Pn- 



e. 



Y,e-^'U,Qa{Q,Pr.-)-k) 

4=1 

(fc) 



fcE 



E^~'^' 



,i=l 



D 



OPTIMAL EXECUTION WITH PRICE IMPACT 23 



Theorem 4.2. V''^'' solves the Hamilton-Jacobi-Bellman 



dip d(p 



mm ^ Pip - Aip,j{p)-^ + -^-p} ^0, ( |3.4 | 



with 7 as in (3.2 1. Therefore, by theorem 3.2 V^^' — Vq 



Proof. First, consid er th e case when there is no impact in the price. Then 7 = 
and by proposition [2t1 l^(°)(y) = xU'^°\p), where C/("'(p) == sup^E[e-''^P^]. 
Since U is viscosity solution of 

niin {j3ip — Aip, (^ — p} = 



Then y(") is solution of (pll) 



Assume now that there is price impact. We know that F'"-* satisfies the 
equation 

min \Bip- Aip,ip- M^^'^iP !• = 0. 



lin (l3ip - Aip, if - 7\/f ("VJ 



Supersolution property: Let y^ ^ O and cp £ (7^(0) such that yQ is a minimizer 
ofFi'''-(^onC'withKO(2/o) ^^{yo)- Hence, pvj°\yo) - Aipiyo) > 0. Now, let 
< C < a^o such that M^°'>ip{yo) = ip{xo - C,aiC,Po)) + Ca{C,Po)- Thus, 

0<V.j'Hyo)^M^°W^!"\yo) 

< VJ°\yo) - v}°\xo - C,a{C,Po)) - Ca{C,Po) 

< f{yo) - f{xo - C*,a(C*,Po)) - Ca{C,Po) 
= ^(2/o)-X^°V(yo). 

Since (p < A^'^V, then ip{yo) = M^^^ifiya). This implies that C = is a 
maximum for ( t-^ ip{xo — C, a(CjPo)) + Co^iC^Po), therefore 



C=o 



da , ^ ,dp, , dip, , /> X >9a,, , 

> ^(C,Po)^(yo) - ^(2/0) + «(Cpo) + Cg^iCpo) 
= -iiPo)-g-{yo) - "^(yo) +Po- 

Subsolution property: Let yo £ O and ip G C^{0) such that j/o is a maximizer 
of l/(°)* — ip on O with V'*-'^-'*(2/o) = ^(2/0)- Without loss of generality we can 
assume that ya is a strict local maximum, that is, there exists S > such that 
j/o is maximum of V^^'* — ip over -65(2/0) <^ C- Let (y„) be a sequence in i?5(yo) 
such that j/„ — ^ j/o £^nd 

hm V('Hy^)^V^'>iyo). 
Recall that V'^'^^ is continuous and is the unique viscosity solution of 

mmhip-Aip,ip-M'^''^ip\=0, (4.1) 
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for all fc > 0. Let yk be a maximum of V^''^ ~ (p over Bs{yo) and let y' be a 
limit point of (y^) as A; ^ 0. For all fc and all n wc have that 



By lemma |4.1[ taking A; — > along the sequence such that yk — > y', we have 
that for all n 

Taking n — >■ oo we obtain that 

and therefore y' = yo since yo is a strict local maximum. Thus, yu — >■ yo as 
fc — >• 0. Let e > , since T^(°)(y„), (/3(y„) — >• (pijjo)^ then for a fix n large enough 
v(2/n) ~ y^'^\yn) < e and for k small enough T^^°''(y„) — V'^^^yn) < £■ Hence 

vivk) - l^(")(yfc) < ^(y„) - l^(''(2/n) < 2e. 

The above shows that y^ is a local maximum of y^*"'^ — ip'^^'^ over Bs{yo) where 
p-('=)(yfc) = Lp(^\yk), (^e^) = ^ - efc and < efc ^ as fc ^ 0. Since F^^) is 
subsolution of ( |4.1[ ), we can consider two cases: 

• There exists a sequence such that (3(p^''^yk) — Aip'-'^^^yk) = (3V'^''^{yk) — 
Aip'-''\yk) < 0. Taking fc — ^ along the sequence we have that 

/3^^°^*(yo) - Aifiiyo) = Myo) - Aifiiyo) < 
and X^'"^ is subsolution. 

• For all k small enough j3Lp''^\yk) — Aip'^^\yk) > 0. This implies that there 
exists < Cfc < 2;fe such that 

XF('=)(yfe) = yW(xfe-a,a(a,Pfc))+Cfe«(a,Pfe)-A:>¥'^'=Hyfc). (4.2) 

Let C,' be a limit point of (Cj.) as fc -^- 0. We claim that C,' — 0: Suppose 
C > 0, then for k small enough such that < C'/2 < Cfe 

MV^^Hyu) > V^'^Hxk - C72,«(C72,Pfc)) + a(C72,Pfc)C72 - k 

> V^^Hxk - Ck, a{Ck,Pk)) + aiC/2,pk)C/2 + a{a,Pk){Ck - C72) - 2fc 
= MV^'^Hyk) + C/2[a{C/2,pk) - a(Cfe,Pfe)] - k. 

where the second inequality follows from ( |2.4[ ) and the definition of V^^K 
Now, since a is strictly decreasing in ( (otherwise condition ( |2.4| cannot 
hold) we can choose k small such that C72[a(C72,Pfc) ~ ct(Cfeil'fe)] > k 
and we get a contradiction. 

Since (^k — > 0, for k small we have that {xk — Cfcj Q!(C/c,Pfe)) G Bs{yo). 



Therefore, from (lO) we have that (f'^'''' [xk - Ck, a{Ck,Pk)) + Cka{Ck,Pk) 



k > •^^^\yk), that is the same as 

•f{xk-C,k,oi{C,k,Pk))+C,ka{C,k,Pk) > ^{xk-Ck,a{Ck,Pk))+Cka{Ck,Pk)-k > (p{yk)- 

(4.3) 
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If we consider, for each k, the map C H' ipixk — C,a{(^,pk)) + Ca{(^,pk), 



(4.3) imphes that there exists Cfc G (0,Cfe) such that 



< ;7^(4,Pfc)^(a;fc-Cfc,a(4,Pfc))--^(a;fc-Cfc,a(Cfe,Pfc))+a(Cfe,P/c)+Cfc^(Cfe,Pfe)- 



Taking fc ^ 



< ~i{po)-g-{yQ) - 'g^(yo)+Po- 



D 



5 Computational Examples 

We are going to present different choices of price processes. Throughout this 
section we wiU consider the price impact function: 

a{x,p) — pe^^^ 
lip) = Ap 

In the following examples, analytical solutions for V do not seem easy to find. 



so we used an implicit numerical scheme following chapter 6 in Kushner and 



Dupuis 1992 . In particular, we used the Gauss-Seidel iteration method for 



approximation in the value space. Additionally, for the impulse control case we 



followed the iterative procedure described in 0ksendal and Sulem 2005 



5.1 Impulse control with positive fixed transaction cost 

Consider the price process following a geometric Brownian motion with fj- < l3 
so that the value function is finite. Figure 1 shows the contour plot of the 
optimal number of shares the investor need to trade. The figure also shows the 
optimal strategy when the investor starts with 5 shares at a price of 2. At time 
0, the investor needs to trade three times until the state variable enters the 
continuation region C (i.e. when the optimal number of shares is 0).When k is 
smaller, the number of trades at time increases and the continuation region 
shrinks. When A: = we obtain the situation described in theorem 13.31 



5.2 Singular control 

Consider the case when the price process follows an arithmetic Brownian motion. 
Then the price dynamics are 

dPt = fidt + adBt - XPtd^t, 

with (T > 0. In this case the value function is always finite, regardless of /i, 
due to the exponential decay of the discount factor. Since is an absorbing 



26 



M. JUNCA 




NjmbernlahareB 



Figure 1: Optimal number of shares with parameters A = 0.5, /i = 2, c = 1, 
/3 = 4 and fc = 0.2. 



boundary for this process the boundary conditions are given by (2.17). Figure 



2(a) shows the value function obtained by the scheme. 

First note that the conditions of Theorem |3.3| are not satisfied, that is 
U{x,p) ^ xp, and therefore T ^ O, as shown in figure [2(b)] Thus, in this 
case the optimal strategy would be to trade very fast in the trading region until 
the state variable hit the free boundary. The figure also shows how the different 
parameters affect the continuation/trade regions. Now, let's see how the change 
in the parameters of the model affect the value function V. Figure [3 (a) | shows 
that the value function is very sensitive to changes in the parameter A for small 
values but not so much for large values. This behavior is common to both pro- 
cesses GBM (described by theorem |3.3[ ) and BM. This means that the bigger 
the investor (i.e. the larger the price impact) the less sensitive to small changes 
in the value of A. Clearly the value function decreases as the impact increases. 

li j3 — 0, the value function would not be finite for any /i > 0, so small values 
of /3 yield a very large value oiV. As /3 increases the effect in V is diminishing. 
Also, the investor has to act greedily and therefore the trade region approaches 
to O and V approaches to W. 

For /i < it is not optimal to wait at all, so y = W, but as fj, increases 
clearly the value function increases in an almost linear fashion. 

The effect of a in the value function is probably the most interesting one. In 
figure [3(d) I we see that it is beneficial for the investor to have some variance in 
the asset but not too much. An explanation for this is that when the variance 
increases it is more likely for the price process to enter the trading region. On 
the other hand, if the variance is too big, the process can hit too fast. Clearly 
the variance of the revenue increases with a, thus as part of future research it 
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(a) Value function in the BM case with parameters A : 
0.5, /^ = 4, o- = 0.5 and /3 = 1. 
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(b) Continuation-trade region in the BM case. The soUd 
hne shows the contour with parameters A = 0.5, fi = 4, 
a = 0.5 and /9 = 1. In the other hnes only the indicated 
parameter has been changed. 



Figure 2: Value function and continuation-trade region in the BM case. 
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(a) Change in V{5, 2) as A varies and jj, = A, (b) Change in V{5, 2) as /3 varies and fi 
a = 0.5 and /3 = 1. 4, ct = 0.5 and A = 0.5. 





(c) Change in V(b, 2) as ^ varies and A = (d) Change in V{'o, 2) as a varies and fi ■ 
0.5, a = 0.5 and 13 = 1. 4, A = 0.5 and 13 = 1. 



Figure 3: Change in the parameters of the model BM. 
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would be interesting to consider the risk aversion of the investor. 

The second example is when the price follows the OrnsteinUhlenbeck process. 
Then the price process becomes 

dPt = a(m - Pt)dt + adBt - XPtd^t, 

with a, a > 0. As in the case of arithmetic Brownian motion, the boundary 



conditions are given by (2.171, since is an absorbing boundary for this process. 
Figured shows the value function and the continuation-trade region. Again, this 
case does not fit within Theorem |3.3[ so the strategy is similar to the BM case. 
Also, the sensitivity of the function to the parameters is similar to the previous 
case. The only parameter that is exclusive to the mean-reverting process is the 
resilience factor a. As we increase a the value function increases (Figure [5 (d)[ ) 
and the continuation region grows (Figure |4(b)[) . 



6 Conclusions 

The main goal of this work was to characterize the value function of the op- 
timal execution strategy in the presence of price impact and fixed transaction 
cost over an infinite horizon. We formulated the problem using two different 
stochastic control settings. In the impulse control formulation we showed that 
the value function is the unique continuous viscosity solution of the Hamilton- 
Jacobi-Bellman equation associated to the problem whenever the transaction 
cost is strictly positive. The second formulation ruled out any transaction cost 
and admitted continuous singular controls only. In this case we also proved 
continuity and uniqueness of the value function under the viscosity framework. 
The next step, part of future research, would be to find the regularity of the 
value function. Numerical results provided in this paper, at least for the second 
formulation, suggest that the function is more than just continuous and that 
its regularity is related with the regularity of the function U defined in Section 
[2j Although any impulse control is a singular control, in general the expected 
revenue obtained when applying the same impulse control in both formulation 
is different. However, the value function may be the same. In fact, we were able 
to show that this is the case for a special type of price impact and provide the 
explicit solution. The question if this is true in general is still unanswered. This 
is particularly challenging since the subsolution property for the HJB equation 



(2.9 1, when there is no transaction cost, has no information at all. Try to find 
answers is part of future work. Another important conclusion is that the HJB 
equation for the regular control formulation, (??) has not enough information 
to characterize the solution. From an economic viewpoint, it would be impor- 
tant to study the effect of the price impact in hedging strategies and how they 
are different to the strategies obtained in classical models, e.g. Delta-hedging 
in Black-Scholes setting. Include utility functions to account for risk aversion 
is another important extension of this work. Finally, the finite time horizon 
natural extension is currently in preparation. 
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Shares 



(a) Value function in the OU case with parameters A 
0.5, a = 4, a = 0.5, m = 5 and /3 = 1. 
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(b) Continuation-trade region in the OU case. The sohd 
line shows the contour with parameters A = 0.5, a = 4, 
(T = 0.5, m = 5 and /3 = 1. In the other lines only the 
indicated parameter has been changed. 



Figure 4: Value function and continuation-trade region in the mean-reverting 
case. 
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(a) Change in V{5, 2) as A varies and m : 
5, o- = 0.5, o = 4 and /3 = 1. 





(b) Change in V{5, 2) as /3 varies and m = (c) Change in V{5, 2) as m varies and a - 
5, cr = 0.5, Q = 4 and A = 0.5. 4, ct = 0.5, /3 = 1 and A = 0.5. 




(d) Change in V{5, 2) as a varies and A = (e) Change in V{5, 2) as cr varies and a 
0.5, o- = 0.5, m = 5 and /3 = 1. 4, A = 0.5, m = 5 and /3 = 1. 



Figure 5: Change in the parameters of the model OU. 
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A Proof of {MVy < MV* 

Let (p he a locally bounded function on O. Let (2/„) be a sequence in O such 
that (j/„) — > yo ^nd 

lim Mifiivn) = (MLpyiyn). 

n— >-oo 

Since ip* is use and F is continuous, for each n > 1 there exists < Cn "£ ^n 
such that 

The sequence (C„) is bounded (since a;„ — t- xq) and therefore converges along a 
subsequence to ^ e [0,xo]. Hence 

(MipYiyo) = lim Mip{yn) 

n— >-oo 

< linisup A^(^*(i/^) 

n— ^oo 

= limsup(/?*(r(y„,Cn)) + Cna(C«,P«) - k 
<(/p*(r(yo,C))+Ca(C,Po)-fc 
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